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ABSTRACT 

Initial  design  study  of  a  low  altitude  satellite  (LASAT)  communication  sys- 
tem is  performed.  The  use  of  many  inexpensive,  low-altitude  satellites  in  random 
orbits  may  be  a  solution  to  vulnerability  of  current  military  communication  sys- 
tems. Statistical  study  of  orbit  characteristics  is  performed  giving  the  mean 
number  of  satellites  in  view  of  a  ground  station  and  the  coverage  density  of  the 
satellite  array.  Waveform  analysis  is  performed  on  a  coded  and  uncoded, 
orthogonal,  noncoherent,  fast-frequency-hopped  M-ary  frequency-shift-keyed 
signal  in  a  Rician  fading  channel  with  optimum  partial  band  jamming.  An  ana- 
lytical expression  for  the  system  probability  of  bit  error  is  obtained  and  numerical 
results  are  generated  for  various  levels  of  fading,  jamming,  and  diversity.  For- 
ward error  correction  coding  is  applied  using  convolutional  codes  and  Reed- 
Solomon  codes. 
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I.     INTRODUCTION 

Current  U.  S.  military  communciations  are  performed  using  a  small  number 
of  highly  sophisticated  and  expensive  satellites  in  geostationary  orbits.  The  loss 
of  only  a  few  of  these  satellites  will  have  a  serious  impact  on  vital  communi- 
cations. Although  current  technology  has  not  yet  produced  an  anti-satellite 
(ASAT)  system  that  can  reach  geostationary  orbits,  such  systems  may  be  devel- 
oped in  the  future.  In  addition,  ground  tracking  and  control  stations  are  suscep- 
tible to  breakdown  and  attack.  To  close  this  vulnerability  window  many  more 
satellites  must  be  fielded  quickly.  The  problem  is  the  enormous  cost  of  building 
and  launching  these  high-altitude  satellites,  which  is  estimated  at  hundreds  of 
millions  of  dollars  each. 

A  possible  solution  to  this  vulnerability  in  times  of  crisis  is  the  use  of  less  so- 
phisticated, low-altitude  satellites  designed  to  be  launched  quickly  in  large  num- 
bers. These  satellites  would  be  light  and  inexpensive  to  build  and  operate. 
Hundreds  of  these  satellites  would  be  placed  in  orbit  to  provide  full  earth  cover- 
age and  their  sheer  numbers  would  overwhelm  any  ASAT  system.  To  achieve 
reasonable  costs,  the  satellites  will  have  minimum  station-keeping  ability,  if  any. 
Initial  system  analysis  has  indicated  that  about  250  satellites  in  random  orbits 
would  be  needed  to  provide  full  earth  coverage  with  inplace  reconstitution 
[Ref.  1].  The  multi-satellite  system  (MSS)  considered  in  Ref.  1  is  primarly  ex- 
pected to  be  used  as  a  strategic  replacement  to  the  current  geostationary  system 
and  is  designed  to  provide  world-wide  communications.  This  configuration 
would  require  that  the  satellites  communicate  between  each  other  via  a  packet 
switching  system.  The  system  considered  in  this  thesis  is  similar  to  the  multi- 
satellite  system  of  Ref.  1  but  would  be  used  for  short  term  tactical  communi- 
cations.  Therefore,  the  satellite  crosslinks  are  unnecessarv. 


II.     BACKGROUND  INFORMATION 

A.  LOW-ALTITUDE  SATELLITES 

The  main  purpose  of  Low-Altitude  Satellites  (LASAT's)  is  to  handle  crisis 
communications  for  military  commanders.  The  satellite  is.  expected  to  operate 
with  encrypted  data  and  digitized  voice  links  to  local  areas,  possibly  between 
hand-held  computers  in  the  field.  One  advantage  of  LASAT  is  the  small  propa- 
gation delay  as  compared  to  250  ms  for  geostationary  satellites.  This  will  facili- 
tate the  speedy  transmission  of  data,  and  many  data  protocols  designed  for 
terrestial  radio  communications  can  be  modified  for  LASAT  use. 

On  the  other  hand.  LASAT  coverage  is  very  small,  and  the  window  of 
communciations  will  last  for  only  minutes  instead  of  hours.  This  requires  a  great 
number  of  satellites  to  provide  continuous  coverage.  It  is  envisioned  that 
LASAT's  will  be  deployed  randomly  in  many  orbits  without  station-keeping  ca- 
pability and  that  communications  will  be  through  omnidirectional  antennas. 
Such  a  random  system  requires  more  satellites  to  provide  continous  coverage  than 
does  one  with  well-defined  orbits  where  station  keeping  is  required.  Therefore. 
the  number  of  satellites  seen  by  a  ground  terminal  is  a  random  variable,  and  a 
ground  station  signal  may  be  received  simultaneously  by  many  satellites.  Also. 
when  used  for  tactical  communications  without  cross-links,  the  same  ground 
signals  relayed  by  many  LASAT's  become  satellite  multipath  signals  at  the  re- 
ceiving ground  terminals  with  near  identical  power  if  classical  transponders  are 
used.  In  addition,  because  low-gain  antennas  will  be  used  by  the  ground  termi- 
nals and  LASAT's.  reflection  multipath  signals  are  also  present.  The  combina- 
tion ol^  these  multipath  effects  result  in  a  fading  channel. 

B.  THE  FADING  CHANNEL 

Generally,  radio  waves  do  not  take  a  direct  path  from  the  transmitter  to  the 
receiver.  Terrain,  buildings,  vegetation,  and  the  ionosphere  cause  multiple  re- 
flections and  refractions  of  the  signal.    Therefore,  a  communications  channel  is. 


in  most  cases,  a  multipath  channel  with  the  signals  arriving  at  the  receiver  from 
a  variety  of  directions  and  with  different  delays.  Furtheremore,  the  channel  is 
very  dynamic  with  the  movement  of  the  ionosphere,  receiver  and  transmitters, 
and  foliage  causing  a  time  variance  in  the  structure  of  the  medium.  As  a  result 
of  such  time  variations,  the  multipath  varies  with  time.  Small  changes  in  path 
length  can  also  result  in  large  phase  changes  in  the  signal  at  the  receiver  antenna. 

Stein  [Rcf.  2]  describes  the  system  in  phasor  terms,  where  the  observed  re- 
ceived phasor  is  a  vector  sum  of  several  phasors  with  each  varing  randomly  and 
individually  over  a  full  2;:  range.  As  a  result,  constructive  and  destructive  inter- 
ference occurs  and  the  received  signals  will  have  different  amplitudes.  The  fading 
channel,  therefore,  is  characterized  by  a  random  time-variant  impulse  response. 
Thus,  the  transmission  of  a  pulse  will  result  in  a  pulse  train  arriving  at  the  re- 
ceiver with  variations  in  amplitude.  If  a  series  of  pulses  are  sent,  the  pulse  train 
as;  dated  with  each  transmission  will  be  different  as  shown  in  Figure  1  on  page 
4  [Ref.  3  :  p.  704].  These  time  variations  arc  unpredictable  to  the  user  and  may 
be  considered  random  which  implies  that  the  received  signal  can  be  modeled  as 
a  random  process. 

In  the  general  case,  the  received  signal  will  have  both  a  direct  (  or  specular) 
component  and  a  number  of  scattered  (or  diffuse)  components.  The  scattered 
components'  arrival  times  and  amplitudes  will  vary  randomly  about  the  mean 
provided  by  the  direct  component.  As  a  result,  the  envelope  of  the  received  signal 
will  be  similar  to  the  envelope  of  a  non-zero  mean  signal.  This  type  of  channel 
is  described  as  Rician.  The  amplitude  of  the  received  signal  has  the  probability 
density  function  (pdf)  [Ref.  4:  p.  105] 

fA{a)  =— rexp^ —  <I0    —    ,       a  >  0.  (2.1) 

cT         I  2o~      J     \  a    J 

where  a2  represents  the  power  of  the  direct  component,  2c1  is  the  expected  value 
of  the  power  of  the  scattered  component,  and  I0(x)  is  the  zeroth-order  modified 
Bessel  function  of  the  first  kind.   For  the  case  where  there  is  no  direct  component 
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figure   1.        Example    of    a    multipatli    channel    on    transmitted    pulses:     From 
Ref.  3:  p.704. 


and  all  the  received  energy  is  via  the  scattered  components,  the  channel  is  de- 
scribed as  Raylcigh  with  the  pdf  of  the  amplitude  given  by 
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The  fading  characteristic  of  the  channel  may  be  rapid  or  slow  dependent 
upon  the  relative  time  and  frequency  scale  of  the  medium  and  the  signal  band- 
width.    The  coherence  bandwidth  is  defined  as  a  measure  of  the  frequency 


dependence  of  the  channel.  If  the  channel  coherence  bandwidth  is  smaller  than 
the  bandwidth  of  the  transmitted  signal,  the  channel  is  frequency-selective  and 
the  signal  is  severely  distorted.  However,  if  the  coherence  bandwidth  is  larger 
than  the  signaling  bandwidth,  the  channel  is  frequency-nonselective,  and  the  sig- 
nal is  undistorted.  A  slowly-fading  channel  is  one  in  which  the  coherence  time  is 
large  and  the  signal  is  unchanged  over  a  signaling  interval.  On  the  other  hand, 
in  a  rapidly-fading  channel,  the  signal  is  distorted  by  changes  in  amplitude  or 
frequency  during  the  signaling  interval. 

Unusual  conditions  in  the  ionosphere  can  result  in  significant  scintillations 
of  earth  satellite  links.  These  effects  are  highly  frequency  dependent.  In  mild 
scintillation  the  channel  is  typically  Rician  but  at  times  of  strong  scintillation 
(deep  fading)  has  been  observed  to  be  Rayleigh.  The  amount  of  scintillation  is 
dependent  upon  satellite  altitude,  time  of  day.  and  latitude,  with  the  largest  ef- 
fects near  the  nighttime  equatorial  zone.  Further  significant  multipath  effects 
that  would  be  involved  in  a  LASAT  system  are  dependent  upon  terrain  and  the 
speed  of  the  satellite  and  terminals.  The  relative  speed  between  the  satellite  and 
the  ground  station  can  be  significant  with  a  low  earth  orbit  satellite  system,  and 
will  increase  fading  effects.   [Ref.  2] 

Diversity  is  one  technique  that  can  be  used  to  overcome  the  effects  of  fading. 
With  diversity  the  same  information  is  provided  over  two  or  more  independently 
fading  channels.  This  lowers  the  probability  that  any  individual  symbol  will  be 
affected  by  deep  fading  on  each  channel.  There  are  a  variety  of  methods  of  di- 
versity combining  and  they  include  the  following  [Ref.  2]: 

•  Spaced  antennas 

Antenna  spacing  of  ten  wavelengths  or  greater  can  result  in  independent 
channels.  Thib  technique,  however,  is  not  well  suited  to  satellite  UHF  com- 
munications because  of  the  long  wavelengths. 

•  Frequency  diversity 

The  use  of  wideband  frequency  hopping  (fast  hopping)  on  fading  chan- 
nels, with  the  same  information  carried  on  different  hops,  tends  to  result  in 
independent  states  o[  fading. 


•  Multipath  diversity 

A  signal  having  a  bandwidth  greater  than  the  coherence  bandwidth  of 
the  channel  will  resolve  multipath  components  and  thus,  provide  the  the  re- 
ceiver with  several  independent  fading  signal  paths.  This  is  an  inherent 
characteristic  of  direct  sequence  spread  spectrum  systems. 

•  Time  diversity 

The  same  data  is  repetitively  transmitted  at  intervals  that  exceed  the 
coherence  time  of  the  channel  which  results  in  independent  fading  condi- 
tions. 

Coding  and  interleaving  are  other  techniques  which  are  similar  to  time  di- 
versity in  providing  redundancy.  This  is  particularly  true  with  binary  linear 
block  codes.  However,  convolutional  codes  and  other  non-binary  linear  codes 
provide  satisfactory  results.  A  more  detailed  coverage  of  fading  channels  is  pre- 
sented in  Ref.  2  and  3. 

C.     FREQUENCY  HOPPING  SPREAD  SPECTRUM 

The  previous  section  discussed  the  effects  of  the  fading  channel  and  indicated 
some  methods  available  to  overcome  the  resulting  degradation  in  system  per- 
formance. Fast  frequency  hopping  spread  spectrum,  with  the  same  information 
carried  on  different  hops,  is  one  form  of  frequency  diversity  that  may  be  used. 
If  sequential  hops  exceed  the  coherence  bandwidth  of  the  channel,  then  fading  for 
each  hop  frequency  is  independent.  The  selection  of  the  frequency  slots  in  each 
signaling  interval  is  made  pseudo-random  according  to  the  output  of  a  pseudo- 
noise  generator. 

Due  to  rapid  signal  phase  changes  which  occur  within  a  fading  channel  and 
the  difficulty  in  building  and  using  coherent  demodulators  for  frequency  shift 
keyed  systems,  most  fast  frequency  hopping  systems  use  noncoherent,  orthogonal 
modulation  schemes.  A  common  data  modulation  for  FH  systems  is  M-ary  fre- 
quency shift  keying  (MFSK),  were  A"  =  log;  M  bits  are  used  to  determine  which 
of  M  frequency  tones  will  be  used  for  carrier  modulation.  The  tones  are  usually 
spaced  far  enough  apart  so  that  the  transmitted  signals  arc  orthogonal.  This 
implies  that  the  frequency  spacing  be  at  least  \jKTb  Hz.  where  Th  is  the  bit  du- 
ration and  Ts  =  KTb  is  the  duration  of  each  symbol. 


When  fast  frequency  hopping  is  used  with  a  M-ary  FSK  system,  each  of  the 
M  symbols  is  subdivided  into  L  chips.  At  each  chip,  the  MFSK  modulator  out- 
put is  hopped  to  a  different  frequency.  Since  the  chip  duration  Tc  is  shorter  than 
the  data  modulator  output  symbol  duration  Ts.  the  minimum  tone  spacing  for 
orthogonal  signals  is  B  =  \jTc  =  LjKTb.  The  total  system  bandwidth  then  be- 
comes IV  =  mMB  =  mLMjKTbi  where  m  is  an  arbitary  integer  that  defines  the 
hopping  bandwidth. 

A  typical  frequency  hopping  system  is  shown  in  Figure  2  on  page  8.  At  the 
receiver,  an  identical  PN  generator  synchronized  with  the  received  signal  removes 
the  pseudo-random  frequency  translation  by  mixing  the  synthesizer  output  with 
the  received  signal  [Ref.  5].  For  an  orthogonal,  noncoherent  system,  demodu- 
lation is  usually  accomplished  through  square-law.  or  envelope  detectors.  The 
detector  output  can  then  make  a  hard  decision  at  each  chip  interval  as  to  which 
of  the  M-channels  contains  the  signal  and  make  an  estimate  on  all  L  chip  deci- 
sions, or  it  can  sum  the  energy  of  L  chips  on  each  channel  and  make  a  soft  deci- 
sion based  on  the  largest  value. 

Although  the  fast  frequency  hopping  system  is  effective  in  a  fading  channel, 
partial  band  jamming  can  seriously  degrade  the  signal  when  it  hops  in  and  out 
of  the  jamming  band.  The  partial  band  jammer  concentrates  its  power  over  a 
small  fraction  ;/  of  the  system  bandwidth.  On  some  hops,  the  signal  will  be  at  a 
frequency  which  is  jammed,  while  on  other  hops,  the  signal  will  be  at  a  frequency 
where  the  only  interference  is  thermal  noise.  A  system  which  can  measure  the 
amount  of  thermal  and  jamming  noise  power  per  hop  can  provide  this  side  in- 
formation to  normalize  the  measurements  obtained  at  the  detector  outputs 
(chapter  IV).  Worst  case  jamming  occurs  when  the  jammer  optimizes  ?/  to  maxi- 
mize the  system  bit  error  probability.  In  system  design  it  is  normally  assumed 
that  the  jammer  has  knowledge  of  the  system  state  and  is  able  to  optimize  ?/  . 
Therefore,  performance  calculations  provide  an  upper  bound  on  error  rates. 

D.     SUMMARY  OF  RESEARCH 

This  thesis  covers  the  communications  analysis  of  LASAT  system  design. 
Included  in  this  studv  are  orbit  characteristics  of  the  svstem.  and  the  waveform 
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Figure  2.      Typical  Frequency  hopped  spread  spectrum  system 


analysis  of  a  proposed  fast  frequency  hopped  M-ary  frequency  shift  keyed  mod- 
ulation scheme. 

1.     Orbit  statistics 

Chapter  III  provides  a  study  of  the  LASAT  orbit  distribution  and  cover- 
age. The  first  section  oi'  the  chapter  looks  into  the  number  of  satellites  in  view 
of  a  ground  station.  An  analytical  approach  was  used  to  derive  the  mean  value 
and  the  standard  deviation  of  the  number  in  view.  A  computer  simulation  veri- 
fied the  results.  The  second  section  of  chapter  III  discusses  the  coverage  of  the 
LASAT  relative  to  a  ground  terminal.  The  mean  coverage  area  and  the  standard 


deviation  in  which  connectivity  can  be  established  with  other  ground  terminals 
was  found.   A  computer  simulation  verified  the  results. 
2.     Waveform  selection 

Waveform  selection  involves  spread  spectrum  techniques,  carrier  modu- 
lation, error  correction  coding,  and  diversity.  Chapter  IV  presents  the  derivation 
of  the  bit  error  probability  for  uncoded  and  coded  fast-frequency-hopping  M-ary 
frequency-shift-keying  schemes  in  a  Rician  fading  channel  with  optimum  partial 
band  jamming.  Numerical  results  for  various  levels  of  diversity  and  fading  are 
provided  for  both  the  coded  and  uncoded  systems.  Convolutional  codes  and 
Reed-Solomon  codes  were  used  to  provide  forward  error  correction  coding. 
Chapter  V  presents  conclusions  and  recommendations  for  further  research. 


III.     LASAT  SYSTEM 

A.     DISTRIBUTION  OF  LASAT 

Since  the  LASAT  communication  system  is  envisioned  as  consisting  of  a  large 
number  of  low  earth  orbiting  satellites  distributed  randomly  throughout  many 
orbits,  it  is  expected  that  the  number  of  satellites  in  range  of  a  ground  station  will 
also  be  random.  The  distribution  of  satellites  within  view  of  a  ground  station  will 
be  developed  here  both  theoretically  and  numerically.  The  ground  station  is  as- 
sumed to  have  an  omnidirectional  antenna  that  will  reach  satellites  with  an  ele- 
vation angle  of  zero  degrees.  A  general  structure  of  L  satellites  randomly 
distributed  within  an  orbit  altitude  range  /z,  to  h2  will  be  assumed. 
1.     Theoretical  development 

The  total  volume  V  enclosed  by  the  orbital  radii  Rx  and  R2  is 

V  =  ±Ti{Rl-Rl)  (3.1) 


where 


R,  =  /?    +  h-> 


(3.2) 


with  Re  =  6378.14  km.  (3443.93  NM),  the  radius  of  the  earth.  The  volume  den- 
sity of  L  satellites  within  this  volume  is 

p  =  L/V     satellites 'km3.  (3.3) 

The  volume  of  the  orbit  portion  seen  by  a  ground  station  can  be  derived  with  the 
help  of  Figure  3  on  page  12.  The  volume  of  the  spherical  sector  vsec2  defined  by 
points  ACBO  with  radius  R:  is 
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where 


»/ 


?.  /»->- 


vsec2  — 


(*x 


r  sin  8ddd(j)dr 


J0     ^0     n) 


—  nR2{l  —  cos  a) 


3 


nRiKRc  +  h2)  -  RJ 


2        2 


(3.4) 


cos  a  =  RJ{Re  +  h2). 


(3.5) 


The  right  circular  cone  defined  by  the  points  ABO  has  a  volume  of 


cone 


fRca- 

f  {R\  -  R2e)Re 


(3.6) 


where  a  is  the  radius  of  the  base  and  a1  =  R;  —  R;  .  The  volume  of  the  segment 
and  base  ABC  is  the  difference  between  the  volumes  o[  the  spherical  sector  and 
the  right  circular  cone. 


v2  -  vsec2  ~  vcone 

=  -j  n{Re  +  h2Y  -  -j  n[{Re  +  h2)2  -  R*~]Re 

=  y  n\2R]h2  +  ARM;  +  2hj  -  {R<  +  2Reh2  +  h2Re 

=  -1-  7ihi[3R2  -3/z?  +  2/??] 

=  4rnh2{3R2-h2). 


-Re)! 


(3.7) 


By  a  similar  method,  the  portion  of  volume  seen  by  the  ground  station  below  h, 
is 
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Figure  3.      Geometry  of  satellite  orbit  distribution. 


v,=   fll^-h^ 


(3.8) 


and  the  subvolumc  v  of  the  orbital  shell  of  thickness  h2  —hi  seen  by  the  ground 
station  is 
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V  =  vs  —  V 


(3.9) 


=  y  nOh\R2  -  hi)  -  -j-  7r(3/?fo  -  h\) 

=  4"  n{3h;R-,  -  3/zfiR,  -  hi  +  h\) 

=  y  k[3A|(/?«  +  h2)  -  3hf{Re  +  hx)  -  h\  +  /zj ) 

=  \n\3Re{hl-h\)  +  2hl-2hll 


Therefore,  on  average,  the  number  of  satellites  seen  by  the  ground  station  will  be 
the  density  of  satellites  times  the  subvolume  v,  or  by  using  (3.3)  and  (3.9) 

nave  =  vp  =  vLfV.  (3.10) 

The  above  result  can  also  be  shown  statistically  with  the  help  of  the  de- 
velopment by  Larson  [Ref.  6:  pp.  143-149].  The  distribution  of  satellites 
throughout  the  volume  V  is  random  and  uniform  and  one  can  assume  that  a 
volume  v  can  be  divided  into  m  non-overlapping  subvolumes  of  size  Av.  where 
Av  =  vjm  is  very  small.  It  is  therefore  reasonable  to  assume  that  there  will  or  will 
not  be  a  satellite  in  each  of  these  small  volumes  and  that  the  probability  that  ex- 
actly one  satellite  will  be  in  each  subvolume  is  pAv  =  pv/m  for  each,  where  p  is 
the  density  of  L  satellites  distributed  throughout  the  total  volume  V.  Further- 
more, it  can  be  assumed  that  the  presence  of  a  satellite  in  each  subvolume  is  in- 
dependent of  the  other  subvolumes.  Thus  the  volume  v  has  been  subdivided  into 
m  repeated  independent  Bernoulli  trials,  each  of  a  volume  Av  =  vjm  with  the 
probability  of  a  satellite  being  present  equal  to  pAv  =  pv/m.  If  n  is  the  number 
of  satellites  in  the  subvolume  v.  then  n  is  approximately  a  binomial  random  vari- 
able with  probabiltity  density  function  fN(n)  where 


A<")  =  (J(4r)(i-Sr)     '    "  =  0A-2 "'         (3.ii 

=  0.       otherwise. 
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If  the  expression  is  taken  to  the  limit  m  -+  co  (thus  the  individual  subvolumes  are 
shrinking  to  zero),  the  exact  probability  function  for  the  Poisson  random  variable 
results.   That  is 


(3.12) 


(PVT      -pv 


Therefore,  the  distribution  of  the  number  of  satellites  within  a  volume  v  can  be 
represented  by  the  volumetric  Poisson  distribution  with  a  parameter  pv.  The 
well-known  mean  n  and  variance  a\  of  the  Poisson  process  is 

n=pv  =  ^-,  (3.13) 

and 

ol  =  f>v  =  ^f.  (3.14) 

2.      Numerical  simulation 

To  generate  a  sample  set  representative  of  the  distribution  over  the  orbital 
volume,  a  uniform  random  sequence  of  L  points  between  0  and  V  is  generated. 
To  find  the  number  of  satellites  seen  by  the  ground  station  for  a  particular  sample 
set  of  L  points  between  0  and  V,  the  number  of  points  between  0  and  v  are 
counted,  where  V  and  v  the  orbital  volume  and  subvolume  defined  by  (3.1)  and 
(3.9).  When  a  large  number  of  sample  sets  are  run  and  the  results  averaged,  the 
mean  number  of  satellites  h  and  the  variance  o;  of  the  samples  can  be  found  from 


M 

A 


and 
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M 


a2 


M 


fa-fl)2. 


(3.16) 


1=1 


where  M  is  the  number  of  sample  sets  and  n,  is  the  number  of  satellites  found 
within  the  subvolume  v  in  the  i  th  sample  set. 

Sample  means  were  found  by  generating  first  2000,  and  then  5000  sample 
sets  of  250  and  150  points  uniformily  distributed  as  described  above.  The  num- 
ber of  satellites  found  in  the  sub-area  were  then  averaged  and  compared  against 
the  mean  found  by  equation  (3.13).  Table  1  shows  the  numerical  results  for 
/?,  =  550  km  (300  NM)  and  h2  =  750  km  (400  NM)  with  2000  and  5000  samples. 
The  numerical  results  compare  favorably  with  the  theoretical  results  obtained  by 
(3.15)  and  (3.16).  Differences  are  due  to  the  pseudo-randomness  of  the  uniform 
number  generator  used. 


Table   1, 


MEAN  AND  VARIANCE  OF  SATELLITE  DISTRIBUTION. 


Theory 

2000  Samples 

5000  Samples 

L 

n 

cl 

n 

<7;, 

a 

c: 

250 

11.54 

11.54 

12.51 

10.43 

12.48 

10.89 

150 

6.92 

6.92 

7.89 

6.78 

7.87 

6.52 

When  an  experiment  described  by  a  Poisson  process  is  repeated  a  great 
number  M  of  times,  then  it  can  be  expected  that  the  number  Mn  of  times  that 
exactly  n  satellites  are  observed  will  be  [Ref.  7:  p.  158] 


Mn*Mfs{n). 


(3.17) 


To  verify  this,  the  number  of  satellites  found  within  the  subvolume  for  the  run 
of  M  =  5000  sample  sets  was  plotted  versus  the  number  of  times  (occurences)  Mn 
that  a  particular  number  of  satellites  were  within  view  (sec  Figure  4  on  page 
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17).    This  curve  is  compared  to  a  Poisson  distribution  with  the  parameter  pv  as 
given  by  (3.13),  and  shown  in  Table  1  (see  Figure  4  on  page  17). 

Table  2  shows  the  maximum  number  of  satellites  /^  seen  by  the  ground 
station  with  a  .999  and  .990  probability.  Since  the  distribution  is  Poisson,  the 
probability  P  that  there  will  be  no  greater  than  n^  satellites  seen  by  the  ground 
station  is 


P  =  Pr[«  <  nmax) 


'Vnax 


(pyf 


-pv 


(3.18) 


n\ 


«=o 


The  probability  P0  that  a  ground  station  sees  no  satellites  overhead  is  (3.18)  with 
'W  =  0.  or 


P  -  P~pv 


(3-19) 


and  is  also  siven  in  Table  2 


Table  2.       NUMBER  OF  SATELLITES  SEEN  BY  GROUND  STATION 


L 

"mac 

P  =  0.999 

'^inax 

P  =  0.99 

n  =  0 

Po 

250 

2° 

20 

0.00001 

150 

16 

14 

0.001 

The  above  analysis  shows  the  expected  number  of  satellites  a  ground  sta- 
tion may  communicate  with.  For  a  transmitting/receiving  pair  of  ground 
stations,  the  number  of  satellites  that  may  be  used  for  communicating  is  depend- 
ent upon  a  joint  pdf  between  the  number  of  satellites,  the  orbital  altitude  and  the 
distance  between  the  stations.  An  analysis  of  the  satellite  coverage  will  be  cov- 
ered in  the  next  section. 
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Figure  4.       Distribution  of  the  number  of  satellites  within  view:     Compared  with 
ideal  Poisson  distribution  with  parameter  pv  =  11.54  . 
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B.     COVERAGE 

1.     Theoretical  development 

Now  that  the  mean  number  of  satellites  seen  by  a  ground  station  has  been 
established,  it  is  desirable  to  develope  a  distribution  for  the  coverage  area  of  the 
satellites.  A  satellite  on  the  horizon  of  the  base  station  has  a  sub-point  that  is  an 
angle  a  away  from  the  base  station  (see  Figure  5  on  page  19),  where 

x  =  cos~\RJR2).  (3.20) 

The  satellite  at  point  S  can  communicate  an  additional  distance  of  arc  length 
Rea  to  a  user  terminal  at  point  UT.  Therefore,  the  base  station  at  point  BS  has 
a  maximum  communication  range  of  arc  length  2<y.Re  out  to  the  user  terminal 
through  the  satellite.  This  gives  a  maximum  coverage  area  defined  by  the  solid 
angle  of  Ay.  or  4  cos~l(RJR2)  steradians,  which  for  an  average  altitude  of 
h  =  (h2  +  /?i)/2  =  650  km  (350  NM)  is  an  area  of 

A  =  R;cos-\RJ(Re  +  h)-\ 
=  17.63  xlO6    km2. 

Figure  6  on  page  20  shows  an  example  of  the  coverage  area  provided  by 
a  number  of  satellites.  The  spherical  caps  of  area  covered  by  the  satellites  are 
shown  as  circles  of  an  equivelent  area.  In  this  diagram  there  are  three  satellites. 
SI,  S2.  and  S3,  uniformly  spaced  over  the  total  view  area  of  the  base  station 
(spherical  cap  o[  a  radius  of  arc  length  xRc  from  the  base  station).  Each  satellite 
has  a  range  of  arc  length  aRe  also.  User  terminals.  Tl.  T2.  T3.  T4.  and  T5.  are 
also  randomly  spaced  over  an  area  out  to  a  range  of  arc  length  2y.Rc  from  the  base 
station.  One  can  see  from  this  Figure  6  that  the  user  terminals  have  differing 
numbers  of  satellites  through  which  they  can  communicate  to  the  base  station. 
In  this  example,  the  user  terminals  can  communicate  as  follows:  T4  through  three 
satellites.  T2  through  two  satellites.  Tl  and  T5  through  one  satellite  each,  and 
T3  has  no  satellite  in  range.  In  general,  as  the  range  between  the  user  terminal 
and   the   base   station   increases,   there   are   fewer  satellites   through   which    to 
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BS  =  Base  station 
UT  =  User    Terminal 
S     =  Satellite 


Figure  5.      Geometry  of  satellite  coverage. 

communicate.    If  more  satellites  were  added  to  the  system,  the  communication 
ablility  would  increase. 

The  coverage  density  for  the  total  number  of  satellites  L  in  the  orbital 
volume  V  needs  to  be  determined.  This  density  represents  the  number  of  satel- 
lites that  can  communicate  with  both  a  user  terminal  and  the  base  station  as  a 
function  of  the  distance  r  between  them  for  a  given  number  n  of  satellites  that  arc 
within  view  of  the  base  station.  The  satellites  are  uniformly  distributed  through 
the  orbital  volume  and  their  different  altitudes  will  result  in  coverage  areas  of 
different  sizes.    This  adds  another  random  variable  to  the  process.    Since  the 
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Figure  6.       Example  of  satellite  communication  range. 

altitude  difference  is  assumed  to  be  small  compared  with  the  coverage  area  of 
each  satellite,  to  simplify  calculations  without  greatly  affecting  the  results,  it  is 
safe  to  assume  that  all  the  satellites  are  at  the  same  altitude  and  that  the  orbits 
are  circular.  The  mean  altitude  h  will  be  used  in  calculations.  The  area  covered 
by  the  satellites  can  be  modeled  as  an  X-Y  cartesian  grid  with  the  base  station 
at  the  origin.  The  positions  of  the  individual  satellites  arc  represented  by  coor- 
dinate points  x,  and  r,,  which  are  uniformily  distributed  from  —  a.Re  to  aRe  along 
the  X  and  Y  axes.   Therefore,  the  pdfs  for  the  random  variables  x  and  y  arc 


fx(x)=\/2a.Re, 
fy{x-)=l/2*Re, 


aRe  <  x  <  aRe 

u.Re<y<aRe. 


(3.22) 


with  a  zero  mean  and  a  second  moment  E[x]  =  E[y]  =  o\  =  o;,  =  {aReY/3. 

For  a  given  number  of  satellites  n,  there  are  n  sets  of  randomly  distributed 
coordinates  describing  the  satellite  positions,  x,j\.i=  1,2, ...  ,n  .  Each  o[  the 
satellites  covers  an  area  defined  by  the  radius  of  arc  length  a.Re  around   the 
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satellite  subpoint.     The  average  distance  from  the  base  station  to  a  satellite 
subpoint  can  be  described  by  two  new  random  variables 


zy 


=  7rZA"< 

f  (3-23) 

• =  IT  Z/'-' 


i=i 


where  zx  and  zy  are  the  average  X  and  Y  axes  coordinates  for  the  satellites.  Since 
n  is  a  Poisson  random  variable,  the  distribution  of  zx  and  zy  will  be  conditioned 
on  n.   The  conditional  probability  density  functions  of  zx  and  z  are  given  by 


(3.24) 


where  the  *  is  the  convolution  function.  The  central  limit  theory  provides  that  if 
n  is  sufficiently  large  (greater  than  five),  the  sum  of  the  random  variables  ap- 
proach a  Gaussian  distribution.  Thus.  zx  and  zy  can  be  approximated  as  inde- 
pendent ^identically  distributed  Gaussian  random  variables  with  zero  mean  and 
a  variance  of  c] 


fz,  \  n(zx  I  n)  =       ,-L—  expj  — f-  J 


x.-2tt(7:  I    2a_-    J 


(3-25) 


where 
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a]  =  Eizl]  =  Eizv] 


=  E 


n      n 


n      n 


t  YLx'x>  =  ~r  YLE[x^] 


1=1  J=  I 


(3.26) 


Since  zx  and  zy  are  independent,  their  joint  density  function  is  given  by 
fzrz>.  I  Mx--y  I «)  =  fzr  I  n(zx  I  nJ/z,  I  «(z>- ' ") 


2tk7: 


T  exP" 


r  2,2 

f  zx  +  ^V 


(3.27) 


2a: 


The  base  station  can  communicate  with  an  area  in  the  shape  of  a  circle 
with  radius  2y.Re  and  the  distribution  can  be  represented  as  a  function  of  the 
distance  r  of  the  user  terminal  to  the  base  station.  To  do  this,  define  a  new  ran- 
dom variable 


\/  zx  +  zy   ' 


(3.28) 


which  is  the  ranee  from  the  base  station  to  the  user  terminal,  and 


r*\Mn) 


^(2x2y)BC 


fzxZy\n^zy\n)dzxdzy 


2-o:  * 


\  -;  +  ^  & 


exp< : —  >dzvdz,.. 

I  2cZ      J 


'X  "">•• 


(3.29) 


The  region  C  is  a  circle  o[  radius  r.  The  above  integral  can  be  easily  solved  using 
polar  coordinates  [Ref.  S:  p.  9].    Let 


zx  =  r  cos  d 
zy  =  r  sin  0 
dzxdzy  =  rdrdO. 


(3.30) 
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then 


^i«H")  = 


i 


inc.  ^o    ^o 
f    ~r 


f         r2    1 
r  exp\ >drdd 

l     2o:  J 


1  —  exp 


I    2(7? 


2^-1 


and 


■W'- 1 ») 


*W'B)) 


r  >0. 


<fr 


TexP   — 

I   2(7. 


r>0. 


(3.31) 


(3.32) 


If  the  substitution  for  a]  from  (3.24)  is  made,  the  result  is  the  exponential  density 


F    -3m    1 


r        i     \      \  3/2  —3/7? 

//?  i  «(r  > /?) = T  expi i"  r 


/->0. 


(3.33) 


The  unconditioned  density  function  fR(r)  can  be  found  by  multipling  by  (3.12) 
and  integrating  with  respect  to  n 


/r{}')  =        fR\n{r\n)fN(n)dn 


i 


5n 


exp- 


-  5  77? 


(pv 


(3.34) 


J0     2(a/y2  [  2(*/y2  J 


e~pV«. 


2.     Numerical  simulation 

A  numerical  simulation  of  the  coverage  problem  was  run  to  verify  the 
above  results.  In  the  model,  the  average  number  of  satellites  that  were  in  view 
of  both  a  user  terminal  and  the  base  station  were  found  for  various  ranges  from 
the  base  station.  This  was  done  by  using  a  sampling  raster  of  points  spaced  45 
degrees  apart  on  concentric  rings  around  the  base  station.  For  each  test  point, 
the  number  of  satellites  that  were  within  its  view  range  was  counted  (a  satellite 
was  within  ranee  if  the  distance  to  it  from  the  base  station  was  less  than  c/.Rr). 
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A  total  was  found  for  all  the  test  points  on  each  ring  and  the  number  was  divided 
by  the  number  of  test  points  on  the  ring  (nine)  to  get  the  average  number  of  sat- 
ellites that  could  be  seen  from  that  range. 

The  sub-satellite  points  were  generated  by  using  a  uniform  random  num- 
ber generator  to  provide  X  and  }'  axis  points  that  were  within  a  radial  distance 
of  one  from  the  base  station  (the  value  of  aRe  was  normalized  to  one).  The  set 
of  test  points  were  generated  in  concentric  rings  about  the  origin  out  to  a  radius 
of  two.  For  each  test  point,  the  number  of  sub-satellite  points  within  a  radial 
distance  of  one  were  counted  (since  all  sub-satellite  points  were  within  a  radial 
distance  of  one  of  the  base  station  it  was  unnecessary  to  test  the  base  station  to 
satellite  distance).  The  total  number  of  satellites  in  view  for  each  test  ring  was 
then  averaged  to  provide  an  average  number  of  satellites  that  can  communicate 
with  both  the  base  station  and  the  user  terminal  for  that  range. 

The  test  results  for  three  different  runs  of  13  satellites  each  are  plotted, 
along  with  the  theoretical  result  from  equation  (3.30).  in  Figure  7  on  page  25. 
The  theoretical  results  were  multiplied  by  the  the  number  of  satellites  in  the  array 
to  normalize  them  with  the  simulation  count.  The  exponential  character  of  the 
numerical  density  is  clearly  seen  and  matches  very  closely  to  the  theoretical  re- 
sults. The  difference  out  near  the  edge  of  the  range  is  due  to  the  approximation 
of  the  Gaussian  function  by  the  central  limit  theory  near  the  tails  of  the  distrib- 
ution. 
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IV.     WAVEFORM  ANALVSIS 

A.     UNCODED  FH/MFSK  IN  A  RICIAN  FADING  CHANNEL  WITH 
PARTIAL  BAND  JAMMING 

The  LASAT  communications  system  under  consideration  must  be  able  to 
overcome  a  variety  of  performance  degrading  phenomena.  Chief  among  these 
are:  Additive  White  Gaussian  Noise  (AWGN),  jamming  and  interference,  and 
fading.  Partial  band  jamming  is  the  optimal  jamming  strategy  most  likely  to  be 
encountered.  Interference  may  be  a  result  of,  among  other  things,  other  com- 
munications systems,  multiple  access  strategies,  or  multiple  transmitters  in  the 
link.  For  the  LASAT  communications  system  the  most  likely  fading  encountered 
will  be  Rician  with  a  direct  line-of-sight  component  and  a  random,  faded  com- 
ponent. 

The  communications  system  waveform  under  consideration  in  this  thesis  will 
be  required  to  overcome  the  above  conditions.  The  system  chosen  for  study  is  a 
fast-frequency-hopped  M-ary  frequency-shift-keyed  (FH/MFSK)  system.  It  is 
assumed  that  the  channel  will  be  jammed  by  an  intentional  jammer  whose  power 
resource  is  AWGN.  Channel  fading  is  assumed  to  be  Rician  fading  with  a  ran- 
dom faded  component  due  to  terrain  reflection. 

Lindsey  [Ref.  9]  studied  a  MFSK  system  with  diversity  in  a  fading  channel 
and  showed  that  diversity  is  effective  in  overcomming  the  effects  of  fading  on  the 
channel.  Lee  et  al.  [Ref.  5,  10]  considered  a  special  case  of  the  Lindsey  system 
with  the  diversity  applied  by  fast  hopping  the  MFSK  signal.  Jamming  and  cod- 
ing were  applied  in  the  analysis  but  fading  was  not.  The  study  done  here  is  a 
continuation  of  the  work  by  Lindsey  and  Lee  et  al..  by  considering  jamming  in  a 
Rician  fading  channel  oi^  a  FH  MFSK  system. 

The  transmitter  receiver  structure  similar  to  Fig.  2  of  Ref.  5  is  used.  A  block 
diagram  of  the  transmitter  is  shown  in  Figure  8  on  page  28.  The  binary  data  is 
provided  at  a  rate  Rb=  ljTb,  where  Tb  is  the  data  bit  length,  to  a  M-ary  source 
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coder  which  provides  M  =  2K  symbols  at  a  rate  Rs  =  RJK  Hz.  There  are  L  hops 
per  symbol.  Each  of  the  symbols  is  then  hopped  at  a  rate  Rh  =  LRS  Hz.  Fre- 
quency hopping  is  employed  over  an  allowed  system  bandwidth  IV  Hz.  The 
carrier  is  hopped  to  a  new  frequency  at  least  once  per  symbol  duration.  Thus,  the 
bandwidth  of  the  transmitted  radio-frequency  pulse  is  B  =  2/t,  where 
t  =  TJL,  Ts  =  symbol  duration  and  the  bandwidth  of  the  M-ary  cluster  will  be 
MB  =  2.U/t. 

The  receiver  is  shown  in  Figure  9  on  page  29.  The  dehopped  signal  r(t)  is  fed 
into  M  channels.  The  output  of  the  square-law  detectors  are  sampled  once  every 
r  =  TJL  seconds  to  produce  the  quantities  {xjk},  i  =  1,2, ...  ,M\  k  =  1.2, ...  ,L.  The 
detector  output  samples  are  multiplied  by  the  weights  \jaj  ,  which  are  outputs 
of  a  noise  only  channel  estimator,  to  produce  {z,k}  and  arc  summed  over  the  index 
k  for  each  i  to  form  the  decision  statistics  {-}.  [Ref.  5] 

If  it  is  considered  that  any  portion  of  the  system  bandwidth  W  has  an 
equally-likely  possibility  of  being  jammed  over  the  jamming  bandwidta  Wfrj, 
where  ?/  is  the  jamming  fraction,  then  a  single  cell  of  bandwidth  B  has  a  proba- 
bility //  of  being  jammed  and  a  probability  1-  //  of  not  being  jammed.  Therefore, 
after  dehopping.  the  received  signal  in  the  A'th  hop  can  be  represented  as 

,  ,      fe(0  +  «*(0 +A(0        with  probability    ?/      } 
{sk(t)  +  /:,(/)  with  probability  1  -  ?/J 

where  nk{t)  andy'/:(r)  are  the  independent  thermal  noise  and  jamming  noise  com- 
ponents, respectively. 

Without  loss  of  generality  it  is  assumed  that  the  Arth  hop  signal  is  present  at 
channel  1  and  is  represented  by 

sk{t)  =  akcos{(Dlt  +  6k).  (4.2) 


and  has  a  transmitted  energy  in  the  Aih  hop  given  by 
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8.      FH/MFSK  transmitter:     Figure  1  from  Ref.  5. 
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/"T 

sk(t) 

X 

ak 

Ek  =  E  =  ±         J^L    \it,       ,/c=l,2,...,L,  (4.3) 


where  it  is  assumed  that  all  transmitted  signals  have  equal  energy  E  .  The  pa- 
rameter ak  in  (4.2)  is  the  Rician  variable  representing  the  signal  amplitude.  The 
signal  sk(t)  has  both  diffuse  and  direct  (specular)  components.  Therefore,  as  dis- 
cussed in  section  2.B  the  coefficients  are  Rician  distributed  and  have  the  density 
function  fAk{ak)  given  by  (2.1) 

ak         {       at  +  c.7.  I     /  a, ,/j.i.  \ 

/AkM = -f  exH  -  -i-^±  M  -¥■  (4-4) 

<7  I  2(7"        J      V      c~      J 

where 


y- 


is  the  modified  Bessel  function  o[  the  first  kind  of  order  a  .  In  (4.4),  ak  is  con- 
sidered to  be  the  strength  of  the  direct  component  of  sk(t)  and  2o2  is  the  mean- 
squared  value  of  the  diffuse  component  in  sk(t).  For  Convenience  define 
y\  =  a|/2a2  as  the  direct-to-diffuse  power  ratio  of  sk{t).  If  y,:  approaches  zero 
(absence  of  direct  components),  the  density  fA(ak)  represents  Rayleigh  fading, 
while  if  y%  approaches  infinity  (presence  o[  direct  components)  the  density  func- 
tion represents  a  Gaussian  channel  with  mean  ak.  [Rcf.  9] 

Thermal  noise  is  assumed  to  be  present  in  each  channel  and  is  additive  white 
Gaussian  noise  (AWGN)  with  zero  mean  and  a  variance  of  ay  =  N()B.  where  \u 
is  the  one-sided  (AWGN)  power  spectral  density  and  B  is  the  channel  bandwidth. 
The  noise  is  assumed  to  be  statistically  independent  and  identically  distributed  in 
each  channel.  For  no  jamming,  the  received  /:th  hop  signal  may  be  represented 
as 
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rk{t)  =  akcos(colt  +  dk)  +  nk(t)i     k=  1.2...  X.  (4.6) 

The  signal  can  be  expanded  into  in-phase  and  quadrature  components  which 
gives 

rk(t)  =  ak{  cos  co{t  cos  6k  —  sin  a^f  sin  Bk)  +  nck(t)  cos  colt  —  nsk{t)  sin  oj{t 

=  [ak  cos  tf/c  +  nck{t)~\  cos  w^  -  [_ak  sin  0^  +  /^.(r)]  sin  (axt  (4.7) 

=  xc/c(f)  cos  a)]/"  —  -v^O  sin  co^t, 

where  nck(t).  and  rc5/t(f)  are  the  independent  noise  quadrature  components  in  all 
channels.    The  output  samples  of  the  square-law  detectors  on  the  /cth  hop  are, 

dropping  the  time  dependence,  for  the  signal  channel 


x\k  -  xck  +  xsk 

=  {ak  cos  0k  +  nck)2  +  {ak  sin  6k  +  nsk)2, 


(4.8) 


and  for  the  other  channels 

**  =  4  +  4,      i  =  2.X...  M.  (4.9) 

When  jamming  occurs,  it  is  considered  to  be  present  in  all  M  channels  and  the 
measurement  channel.  The  jammer  is  assumed  to  spread  its  total  available  power 
J  uniformly  over  the  fraction  ;/  of  the  bandwidth  IT.  0  < //  <  1.  The  average 
jamming  spectral  density  is  Nj=  J\W  Watts  Hz.  The  jamming  spectral  density 
over  the  actual  jammed  fraction  of  the  band  is  Njfrj  Watts.  The  combination  of 
jamming  and  thermal  noise  on  the  kih  hop  produces  the  detector  output  samples 
for  the  signal  channel 

xlk  =  {ak  sin  Bk  +  nck  +jckf  +  (ak  cos  dk  +  nsk  +jskf,  (4.10) 

and  for  the  other  channels 

xu<  =  (nck  +Jckf  +  (>hk  +Lif<       i  =  2,3, V/,  (4.1 1 ) 
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where  jck,  andyjft,  k—  1,2. ...  ,L  are  the  independent  jamming  noise  quadrature 
components  in  all  channels  and  are  assumed  to  be  additive  Gaussian  noise  with 
variance  defined  by  a)  =  NjBjr\  . 

Following  the  derivation  by  Lee  et  al.  [Ref.  5],  the  error  probability  ex- 
pression can  be  derived.  If  equally  likely  M-ary  symbols  are  assumed,  the  ex- 
pression for  the  symbol  error  probability,  given  signal  of  channel  1  is  sent,  is 

L 

ps(e)  =  J]PjPs{e  I  /  hops  jammed)  (4.12) 

/=o 

where 

P^(^\\\-n)L-1  (4.13) 

is  the  probability  that  /  out  of  L  channels  are  jammed,  and 

Ps(e\f)  =  l  -Ps{c\t)  (4.14) 

is  the  probability  that  an  error  was  made  given  /  channels  are  jammed.  P5(c\f)  is 
the  probability  of  a  correct  decision  given  that  signal  of  channel  1  is  sent  and  / 
hops  are  jammed. 

Whalen  [Ref.  4:  sec  4.7.  4.8.]  shows  that  the  square-law  detector  output  has 
a  probability  density  function  which  is  noncentral  chi-squared:  thus,  for  a  given 
amplitude  ak, 


fx^lk  I  **)  =  — T  CXp, -T—      If)     2 

2ck  I  2<7 k       J     l         trA 

where  c\  =  >XB  =  2/V/t  =  2iNLjTs.   The  total  noise  density  N  is  given  by 


(4.15) 


f.\u  with  probability  1  —  ?/l 

A  =  <  >  (4.16) 

I A  /'=  (J^o  +  -V/'/)     w^h  probabilty  ?/.       J 
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The  output  of  the  detector  is  then  normalized  by  dividing  by  the  output  of  the 
channel  estimator  which  results  in  the  variable  zik  =  x]kja}.,  i=  1,2, ...  ,M,  which 
has  a  conditional  pdf.  for  the  signal  channel. 


1 


/z1Jk(-i*K-)  =  yexP- 


-u- 
? 


} 


2*1 


(      ^1/2 
\z\k)     ak 

Ob 


(4.17) 


The  unconditional  probability  density  function  for  zlk  can  be  found  by  integrating 
(4.17)  over  all  ak 


fzj-ik)  =      fz  | Ak{zik  I  aj^A(a^dah 


(4.18) 


or  (see  Appendix). 


i—eJ-="  +  2"*\J    Ul 


/z"(ru)"  2(i-t-/])  expj  ~~W^1 


\k2PI< 


(i+py 


(4.19) 


where  p  =  2o:R  is  the  output  signal-to-noise  ratio  produced  by  the  diffuse  com- 
ponent, and  pk  =  ccj.R  is  the  output  signal-to-noise  ratio  o^  the  specular  (direct) 
component  of  the  k  th  hop  signal  in  channel  1.   The  ratio  R  is  defined  as 


R  = 


_E_ 
N 


hop  energy  of  transmitted  signal 
total  noise  density 


4.V       2B.X 


2a\ 


(4.20) 


The  decision  statistics  at  the  output  of  the  summers  are 


i. 

.  ._  V- 

-l-  IJ-ite 


i=  1.2....  ,M- 


(4.21) 


The  probability  density  function  for  z,  may  be  found  using  characteristic  func- 
tions.   With  the  use  of  the  tables  oi^  Fourier  transforms  by  Campbell  and  Foster 


33 


(pair  655.1)  [Ref.  11:  p.  79]  the  characteristic  function  offz(zlk)  in  (4.19)  is  given 
by 

^)=itrr«^^}7TFexp{ii(^r}-   (4-22> 

where  p  is  the  complex  variable  and 


V^="T(l+/0and   ft 


4(1  +/?)• 


(4.23) 


Since  the  z,/('s  are  independent  and  identically  distributed,  the  characteristic 
function  of  the  decision  statistic  z,  can  be  found  from 


A-l 


Combining  (4.22)  and  (4.24)  yields 


(4.24) 


#z,(/>)  = 


[2(1  +  /I)]L 


exp- 


L 

-Ylo 
r        L2Pk 

\         k=] 


{   2(1+/?) 


1 


fa  +  lW 


cxp 


f 


"l  '■ 


1 


(4.25) 


From  Ref.  1 1:  pair  650.  p.  77,  the  probability  density  function  for  z,  is  obtained 
as  follows: 


L-\ 


fZl(=i) 


f         Z,   -r  2/);    ) 

exp< 


2{l+P)\2pIJ  'l{      2(1+/?) 


•'l-i 


--,2/v 


lV  (l  +  /?r  J 


(4.26) 
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Z]  >  0  .   The  noncentral  parameter  of  (4.26)  is 


Pi 


-Z»- 


fe=i 


L  L 

)^  =  2<72*\yf, 

k=\  k=l 


(4.27) 


where   yjj  =  c/.:J2u2  is  the  direct-to-diffuse  power  ratio  and 


£         Es        KEb 


N       XL        NL 


(4.28) 


is  the  energy  per  hop  to  noise  density  ratio,  with  Es  as  the  symbol  energy,  and  Eb 
as  the  bit  energy.   Then,  if  it  is  assumed  that  /  out  of  L  channels  are  jammed. 


Pi  =  2(7' 


E_ 

A 


y]< 


k=\ 


2c- 


l     KEb        {L-I)    KEb 
L    LNT  +       L        LN0 


2(7 


k=\ 

L   ' 


'Eb/        ,    Eb 

l^-  +  (L-[)—- 

At  ;>n 


(4.29) 


when 


J_ 
L 


y'k 


k  =  \ 


(4.30) 


is  the  mean  direct-to-diffuse  power  ratio  with  A',,  and  NT  are  from  (4.16). 

With  the  use  of  the  definition  of  the  Bessel  function  (4.5).  equation  (4.26)  can 
be  expressed  as 
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fzfa) 


,-pll(l+p) 


c__  ^wi+/y      f    -Zl    1/     ,t     x*" 

2(1+/))    £jp\r(p  +  L)       l  {  2(1+/?)  JV  2(1+/?)  / 


As  in  the  non-fading  case  with  a  fixed  amplitude,  the  pdf  for  the  noise  only 
channels,  z.,  i  =  2.3, ...  ,M,  is  chi-squared  with  2L  degrees  o£  freedom  and  may  be 
obtained  from  (4.32)  with  p,  =  0.    zx  =  z,:,  and  /J  =  0 


(z//2) 


L-l 


/z;fo)  =     ;r(L)     exp(-z,/2).       z,>0,     /  =  2.3,.....U. 


(4.32) 


Equation  (4.14)  can  then  be  found  from 


Ps(e\f)  =  l-       fzXzi 


-.M-l 


fz^dz 


d: 


(4.33) 


Substituting  the  appropriate  distributions  (4.31)  and  (4.32)  into  (4.33)  above 
yields 


Ps(e\f)=  1  " 


pi 

C       l+B 


(p,/2(l  +/?))'      .    f       -r,       )/ 


exp< 


2(1  +/*)  Zj    ^r(p  +  L)       'U  2(1+/))  JV  2(1 


+  « 


2f(L) 


L-l 
-^  e~2'/2<fz 


^0 


A/-1 


«fe. 


p+L-1 


(4.34) 


The  solution  to  the  intesral  inside  the  brackets  is 


-,.V/-1 


L-l 


/z(# 


—  f 


2V(*i/2)' 


,   -,/W-l 


/•: 


r=0 


(4.35) 


which  can  be  expanded  to 
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M-\ 


m=U 


M 


m 


K-iy 


-z,/2V(m/2)' 


r! 


;-=0 


(4.36) 


Af— 1 


»j=0 


'M 


/7Z 


l(- 


mL—m 

hm  -Wi/2    \      C,.{m.L) 
I)  e       " 


r! 


r=0 


9 


(4.37) 


where 


with 


ii" 


Cr(/;z?L)  =  -jr 


[(m+  1)«  -  r]Cr_B(m,L) 


«=i 


(4.38) 


p-        for       0  <  r  <  L  -  1 
[L  -  1      for      r  >  L  -  1 


(4.39) 


Equation  (4.38)  is  a  modified  version  of  the  J.  P.  Miller  formula  [Ref.  12  :  p.  42] 
for  the  coefficients  of  a  composition  of  a  formal  power  series  and  a  nonunit. 
Composition  is  the  substitution  of  one  power  series  into  another.  The  Miller 
formula  applies  when  one  series  is  a  binomial  series.  Substituting  (4.37)  into 
(4.34)  gives 
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Ps(e\f)=l- 


e 

Pi 

y 

1  (p//i 

+  df  i 

1 

+  P 

L 

iP'-f{p 

+  L)  _ 

p=0 
■V— 1  mL—m 

zr:,>-"-x^a 

fM=0 

(1+/?)' 


r=0 

-       r       [m(l+W+l]z, 
exp^  - 


o 


(4.40) 


[  2(1+/?)  JV2(1+/?) 


/7+L+r-l         -. 

dz-i    . 


It  can  be  shown  that  [Ref.  13:  p.  438] 


M-\ 


'M 


mL—m 


m 


CJm.L)        , 


m=0 


r=0 


mL—m 


'  (M  -  1 


m=\ 


m-n    \      CJnhL)        , 


r: 


r=U 


(4.41) 


Solving  the  integral  and  and  using  (4.41),  yields 


M-\ 


_  vY^-i\ (-i 


VW?+1 


m=P     ™     /  (I  +W  +  /72/; 

mL—m 


\     'Pi 

exP1  T— 77 


^  —  11+/? 


C>?.L) 
r! 


r=0 


+  j» 


/>=0 


1  +  m  +  mfi 

(p  +  r  +  L-  1)' 
p!r(/?  +  L)    (l  +/j)/,(l  +m  +  mpy 


W 


(4.42) 


With  some  rearanging  (4.42)  can  be  written  as 


3S 


A/-1 


[M-\\        (-1) 


m+l 


\\    m    /  (i  +  m  +  /7ifl 

mL  —  m 


rexPu  +  /J 


^     Cr(m,L) 


r=0 


1  +  /? 


I  +  m.+  /»/? 
(p/ 


(r  +  L-1)! 
r!(L-l)! 


(4.43) 


>=o 

The  portion  in  the  brackets  is  equal  to  the  confluent  hypergeometric  series 
,F,(/-  +  L:  L:  w)  [Ref.  12:  Vol  II,  p.  260,  Ref.  14:  p.    1250]  and  is  given  byi 


^(r  +  L:  L;  w)  =  -7— =  e^F^  -r,  L\  -  w) 

I  Y  T  ii   —    1 


(4.47) 


when 


w  = 


Pi 


1  +  fi(\  +  m  +  m(3) 


(4.48) 


1  Details,  on  the  derivation  of  equation  (4.47). 

00  oc 

T^  (L  +  \)p    ^p       V^  f (L  +  r  +  />)       f(L)       / 


where 


F,(r  +  L;L:iv) 


11  l 


'        (L)p       p\        /_j     V(L  +  r)       Y(L  +  P)   p\ 

P=o  p=° 


(L  +  r  +  p-W.     /        (L-  l)Kr!)  _£ 

(L  +/>-!)!       (p!)    (L+r-  1)!  (r!)  "   r! 


^\w)jfr+L- 


p=0 


F,(a:b:w)  = 


i1  1 


T(a  +  p) 


p=q 


(4.44) 


(4.45) 


(4.46) 
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L;,"(u-)  is  the  generalized  Laguerre  polynomial,  which  in  this  case  equals 


=  —^-iFi(-r;L;-w). 


(4.49) 


Substituting  (4.47)  into  (4.43)  gives  the  desired  result 


M-\ 


Ps(e\f)=  T 


'A/ 


(-D 


w?+l 


OT=1 


f  ~"*P/  ] 

(1  +  m  +  mP)L         i  1  +  m  +  mP  j 


exp 


1+jS 


mL—m 

;-=0 
r 

"V 

Zj/?!(1  +/?f(l  +  m  +  m$)p 

p=0 


1  -f  /?Z  +  /??/? 


>  +  L  -  1 


(4.50) 


The  result  from  Lindsey  [(41).  Ref.  9]  for  the  symbol  error  can  be  obtained 
from  (4.50)  by  letting  ;/  =  0.  the  no-jamming  case,  and  using  the  hypergeometric 
series  definition  (4.47) 


.u-i 


>m+l 


Pf(c)  =  /  rexn  t: — i — 7 1 

Lj\    ™     /  (1  +m  +  m[J)L         I  l  +  m  +  m0  j 

m=  1 

mL— m 


Cr(m.L)    T{L  +  r) 


r={) 


V(L) 


+  P 


1  +  m  +  /;?/; 


(4.51) 


J 


■1^1 


•r:  L: 


Po 


(1  +/?)(!  +m  +  mP)  _ 


where  p0  is  (4.30)  with  /=  0. 

Bit  error  probability  from  symbol  error  for  an  M-ary  orthogonal  system 
(MFSK)  is  given  by 


41) 


Pb(c)  =  2(UU_1}  P5{e)  (4.52) 

where  the  probability  of  symbol  error  is  given  by  equation  (4.12).    Using  (4.13), 
(4.50).  and  (4.52)  in  (4.12)  gives  the  probability  of  bit  error  for  the  receiver 


2{M-\)  /_j\l)  ^\    m     )  {i+m  +  mp}L 


/=o 

mL—m 

exp'  —  >     >    C.,{nuL)\  

*\  1  +  m  +  m\]  j  l_j     A         \_  1  +  m  +  mfj 

V  iPif  (P  +  L—\ 


-,r 

(4.53) 


/_j  p\{\  +  (if{\  +  m  +  mfif  V     p-r     J 
P=o 

For  the  non-fading  case,  2u2  =  0  and  <xk  —  1  for  all  k,  and  if  the  definition  for  the 
Laguerre  polynomial  (4.49)  is  used  the  result  is  equation  (16)  from  Ref.  5 

^  \     ™     /   (/w  +  1  T         w"  +  1       7 

ml-m  (4.54) 

V     C("?^)    L(L-1)/  Pi 


/  j   (/;/+  1)  V       m  +  1 

r=(J 

B.     CODED  FH/MFSK  IN  A  RICIAN  FADING  CHANNEL  WITH 
PARTIAL  BAND  JAMMING 

The  diversity  studied  in  the  last  section  is  an  effective  means  of  overcoming 
the  detrimental  effects  of  fading.  Frequency  diversity,  however,  is  a  trivial  form 
of  repetitive  coding  and  the  decision  statistic  of  the  receiver  is  soft  decision  de- 
coding of  the  repetitive  code.  More  complex  types  of  codes  may  also  be  used  to 
provide  diversity  which  is  much  more  efficient  than  frequency  diversity  [Ref.  3: 
p.  753]. 
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Forward  Error  Correction  (FEC)  codes  can  be  used  to  provide  improved 
performance  with  fading  signals  without  having  to  use  large  amounts  of  diversity 
which  cause  increases  in  the  recombining  losses  caused  by  the  summing  oper- 
ations. To  apply  coding,  the  circuit  of  Figure  8  on  page  28  is  modified  with  the 
addition  of  a  binary  to  Q-ary  encoder  and  a  Q-ary  to  binary  decoder  as  shown 
in  Figure  10  on  page  43.  Information  bits  from  a  binary  source  at  a  rate  Rb  are 
mapped  to  Q-ary  symbols  where  O  =  2q  at  a  rate  Rq  =  Rjq  and  applied  to  the 
coder  input.  The  forward  error  correction  coder  outputs  n  Q-ary  symbols  for 
every  k  Q-ary  input  symbols.  This  gives  a  code  rate  r  =  kjn  and  the  coder  output 
rate  is  Rc  =  nRJk  =  nR-Jkq.  The  channel  modulation  is  FH/MFSK  with 
M  =  2K.  The  modulator  input  rate  is  Rd  =  qRJK  =  nRJkK.  The  M-ary  words 
are  applied  to  a  MFSK  modulator  which  selects  one  of  M  baseband  frequencies, 
fu  i  =  1,2, ...  ,M.  The  modulator  output  symbols  are  then  hopped  at  a  rate  of  L 
hops  per  symbol,  or  Rh  =  LRd  =  nLRJkK.  The  bandwidth  of  the  transmitted 
pulse  is  B  =  2/t,  where  -  =  TJL  ,  Td  =  symbol  rate.  The  total  system  bandwidth 
is  made-up  of  A^  possible  frequencies  spaced  B  Hz  apart  so  that  W  —  NhB  Hz. 
[Ref.  10] 

The  receiver  block  diagram  is  also  shown  in  Figure  10  on  page  43.  The 
demodulator  operation  is  the  same  as  described  in  section  4. A.  After  a  decision 
is  made  as  to  which  channel  has  the  largest  output,  the  M-ary  words  are  con- 
verted into  Q-ary  symbols  and  sent  to  the  decoder  input.  The  decoder  makes  a 
hard  decision  as  to  what  Q-ary  symbol  was  sent.  A  Q-ary  to  binary  converter 
then  reconstructs  the  binary  series  {X}  . 

If  the  transmitted  energy  per  information  bit  is  defined  by  (4.28),  then  the 
energy  per  coded  M-ary  word  is 


kKEh 


and  the  energy  per  hop  is 


E;       kKEh 

E"  =  —  =  -IT-  (4J6) 
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Figure  10.      FH/MFSK  system  with  coding:     Figure  1  from  Rcf.  10. 
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The  probability  of  word  error  at  the  output  of  the  decision  maker  is  the  same  as 
without  coding,  namely 


M-\ 


m=\ 


M  -  1 


(-D 


m+\ 


m    J  (l  +  m  +  mfi)L 


\        ~mPi 

exp< - 

I  1  +  m  +  mp 

V  (P/ 


mL—m 


^  Cr{m,L) 


r=0 


+  P 


—,r 


1  +  m  +  mp 


P=o 


'p  +  L-1' 


p\{\  +  p'f{\  +m  +  mfif  I     p-r 


(4.57) 


The  noncentral  parameter  pt  for  the  rate  =  k  n  coded  waveform  with  /  out  of  L 
channels  jammed  is 


Pi  = 


kK 
nL 


I-    -  +  (L-f) 


N 


Nt 


(4.58) 


where  NT  =  A"0  +  Njft]  as  before. 
1.     Convolutional  codes 

The  convolution  codes  considered  here  are  binary  codes  with  constraint 
length  7,  and  a  rates  of  1  2  and  1/3.  The  Viterbi  decoding  algorithm  is  assumed 
and  the  information  bit  error  rate  for  the  best  rate,  constraint  length  7  code  can 
be  bounded  by  the  Chernoff  bound.  A  best  rate  convolutional  code  is  one  that 
gives  the  best  performance  from  a  Viterbi  decoder.  [Ref.  15]  For  best  rate  1/2 
codes,  the  bound  is 


^<-^-(36Z)10-r211Z)12+  1404Z)I4+  11633Z>16  +  -  ),  (4.59) 


and  for  rate  1  3 


Pb  <  4"  (7Z>!5  +  SZ)16  +  22D17  +  44Z)18  +  -  ). 


(4.60) 
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where 

D  =  2JPC(1  -  Pc) .  (4.61) 

In  (4.61)  Pc  is  the  forward  transition  error  of  a  binary  symmetric  channel,  which 
for  a  binarv  code  and  a  M-arv-to-Binarv  decoding  is 


ia 


M/2 
M 


Pc  =  -rr^T  ps,  (4-62) 


where  Ps  is  given  by  (4.57).  The  higher  order  terms  for  D  in  (4.59)  and  (4.60)  are 
neglected  due  to  the  usually  small  size  of  Pc 
2.     Reed-Solomon  codes 

Reed-Solomon  codes  were  found  to  be  particarly  well  suited  for  usage  on 
fading  channels,  more  so  than  convolutional  codes  because  Viterbi  decoders  are 
not  well-suited  to  the  bursty  errors  associated  with  fading  [Ref.  16].  The  Reed- 
Solomon  codes  are  non-binary  linear  block  codes  in  which  the  elements  of  the 
code  word  are  selected  from  O  symbols.  Usually  Q  =  2q  so  that  k  information  bits 
are  mapped  into  one  0  symbol.  For  q  >  K,  where  K  =  log2M,  qjK  M-ary  words 
are  mapped  to  the  Q-ary  symbols  and  the  symbol  error  at  the  output  of  the  de- 
coder is 

Pc=\-(\-Ps)q!K.  (4.63) 

In  practice.  q/K  is  usually  an  integer. 

For  an  (k,n)  block  code,  the  decoder  accepts  an  /z-tuple  of  Q-ary  coded 
symbols  and  outputs  a  /c-tuple  of  decoded  information  symbols  which  are  further 
broken  down  into  information  bits  by  the  Q-ary-to-binary  decoder.  At  the  de- 
coder output  the  probability  of  error  for  the  Q-ary  symbols  is 


pq  ^  i  V  (';k(  i  -  Pcr' + 1  i  ("k  i  -  pr\    (^ 


i=i+\ 
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where  d  is  the  minimum  distance  between  code  words,  which  for  Reed-Solomon 
codes  is  d  =  n  —  k  +  1 .  and 


d-  1 


(4.65) 


is  the  maximum  number  of  correctable  symbol  errors  in  a  code  word.  [Ref.    10  ] 
The  notation  [x]  denotes  the  integer  less  than  or  equal  to  x.    The  probablity  of 

bit  error  is  then  given  by 

Oil 
pb  =  -^ZTPr  (4-66) 

C.      NUMERICAL  RESULTS 

in  both  the  coded  and  uncoded  cases,  the  generation  of  the  probability  of  bit 
error  involved  a  computer  simulation  of  equation  (4.53).  In  applying  optimum 
partial  band  jamming,  the  jamming  fraction  ;/  which  maximized  the  probability 
of  bit  error  had  to  be  found.  This  was  accomplished  by  performing  a  numerical 
search  for  ;/  at  each  value  of  jamming  power  to  noise  power  density  ratio  EJNj 
for  every  value  of  nu  L,  EJNQ,  and  ~:. 
1.     Uncoded  performance 

Uncoded  performance  for  various  levels  of  fading  are  shown  in  Figures 
11  to  14  for  EJN0  =  10  dB.  and  M  =  4.  Also  plotted  is  the  ideal  performance  for 
no  fading  or  jamming  with  the  abscissa  as  Eb/N0  vice  EJNj.  In  comparing  these 
curves  for  each  value  of  ?/  one  can  see  the  adverse  effects  fading  has  on  the  system 
performance.  It  can  also  be  noted  how  diversity  improves  performance  signif- 
icantly for  small  value  of  y2  (Rayliegh  fading),  but  the  diversity  also  hinders  per- 
formance when  there  is  a  strong  direct  component  in  the  received  signal  (high 
~:).  At  ~:  =  15  the  performance  of  all  levels  of  diversity  are  about  the  same.  For 
a  signal-to-noise  ratio  of  15  dB.  a  similar  result  is  seen  in  Figures  15  and  16.  For 
the  higher  Eb/N0  .  however,  the  adverse  effects  of  diversity  are  not  realized  in  that 
the  higher  EJNQ  offsets  the  recombining  losses.  Figures  17  and  18  show  a  com- 
parison of  performance  for  various  values  o[  M  when  L  =  2  and  EJN0  =  10  dB. 
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for  both  the  Rayleigh  fading  case  and  a  Rician  fading  case  respectively.    From 
these  curves  it  is  apparent  that  an  increased  symbol  size  gives  better  performance. 
More  improvement  from  diversity  is  seen  in  the  Rician  channel. 
2.     Coded  performance 

The  first  convolutional  code  considered  was  a  best  rate  1/2,  constraint 
length-7  code  with  the  probability  of  bit  error  bounded  by  (4.59).  Plots  of  the 
error  rates  for  this  code  are  shown  in  Figures  19,  20,  21,  and  22  for  various  levels 
of  fading  with  M  =  4  and  EJN0  =  15  dB  in  worst  case  partial  band  jamming.  In 
the  non-fading  and  Rician  fading  with  a  strong  direct  component  cases  (y2  =  15) 
(sec  Figures  19  and  20),  the  optimum  level  of  diversity  is  L  =  2  hops  per  symbol. 
However,  as  the  fading  becomes  stronger  (see  Figures  21  and  22),  the  perform- 
ance increases  when  there  are  a  greater  number  of  hops  per  symbol.  This  is  as 
expected,  if  these  plots  are  compared  with  the  results  for  the  uncoded  system  in 
which  the  higher  levels  of  diversity  increased  performance  as  fading  increased. 
Furthermore,  a  comparision  of  Figures  20  and  22  with  Figures  15  and  16  shows 
that  coding  does  indeed  decrease  error  rates  as  expected.  For  binary-FSK  mod- 
ulation (M  =  2).  the  same  general  pattern  is  evident  (see  Figures  23  and  24)  where 
diversity  improves  performance  as  fading  increases.  Also,  a  comparison  of  the 
BFSK  results  with  the  MFSK  results  indicate  that  better  performance  is  gained 
with  higher  values  of  M  . 

The  best  rate  1  3  constraint  length-7  code  results  are  shown  in  Figures 
25  and  26  for  Rayleigh  fading  and  Rician  fading  when  M  =  2  and  EJN0  =  15  dB. 
Again  one  can  see  that  the  greater  number  of  hops  per  symbol  provide  increased 
performance  in  deep  fading.  With  the  rate  1/3  code,  however,  slow  hopping 
(L=l)  provides  consistently  poor  performance,  even  with  no  fading.  A  (15.8) 
Reed-Solomon  code  was  used  to  generate  the  curves  in  Figures  27  and  28  for 
M  =  4  and  EbjN0=  15  dB  in  Rician  and  Rayleigh  fading  channels  respectively. 
The  probability  of  bit  error  for  this  code  is  given  by  (4.63)  through  (4.66)  with 
d=  8  and  t  =  3.  Therefore,  the  (15.8)  code  has  an  error  correction  capability  of 
3  bits. 
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Figure  1 1.  Performance  of  receiver  for  no  fading:  y2  =  oo  and  M  =  4  under  opti- 
mum partial-band  jamming  conditions  with  number  of  hops-pcr-  sym- 
bol (L)  as  a  parameter  when  EJi\0  =  10  dB. 
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Figure  12.  Performance  of  receiver  for  Rician  fading:  y2  =  15  and  M  =  4  under 
optimum  partial-band  jamming  conditions  with  number  of  hops-pcr- 
symbol  (L)  as  a  parameter  when  EJN0  =  10  dB. 
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Figure  13.  Performance  of  receiver  for  Rician  fading:  y2  =  5  and  M=4  under 
optimum  partial-band  jamming  conditions  with  number  of  hops-pcr- 
symbol  (L)  as  a  parameter  when  EJi\0  =  10  dB. 
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Figure  14.  Performance  of  receiver  for  Rayleigh  fading:  y2  =  0  and  M  =  4  under 
optimum  partial-band  jamming  conditions  with  number  of  hops-pcr- 
symbol  (L)  as  a  parameter  when  EbjN0  =  10  dB. 
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Figure  15.       Performance  of  receiver  for  Rician  fading:     y2  =  15  and  M  —  4  under 

optimum  partial-band  jamming  conditions  with  number  of  hops-pcr- 
symbol  (L)  as  a  parameter  when  EbIN0  =  15  dB. 
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Figure  16.  Performance  of  receiver  for  Rayleigh  fading:  y2  =  0  and  M  =  4  under 
optimum  partial-band  jamming  conditions  with  number  o[  hops-pcr- 
symbol  (L)  as  a  parameter  when  EJi\  =  15  dl3. 
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Figure  17.  Performance  of  receiver  for  Rician  fading:  y2  =  15  and  L  =  2  under 
optimum  partial-band  jamming  conditions  with  number  of  hops-pcr- 
symbol  (M)  as  a  parameter  when  EJN0  =  10  dB. 
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Tigiire  18.      Performance  of  receiver  for  Rayleigh  fading:     y2  =  0  and  L  =  2  under 

optimum  partial-band  jamming  conditions  with  number  of  hops-pcr- 
symbol  (M)  as  a  parameter  when  EJN0  =  10  dB. 
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All  three  codes  discussed  are  plotted  together  in  Figures  29  to  32  for  two 
levels  of  diversity.  L  —  2  (solid  line),  and  L  =  3  (dashed  line),  when  M  =  4  and 
Eb/NQ  =  15  dB.  From  the  curves  one  can  see  that  the  best  rate  1/2  convolutional 
code  consistently  provides  better  detection  capabilities  than  the  other  two.  In 
addition,  the  rate  1/3  convolutional  code  is  slightly  better  than  the  Reed-Solomon 
(15,8)  code  for  signals  with  strong  direct  components,  and  much  better  for  signals 
in  deeper  fading.  Here  again,  the  benifits  of  diversity  can  be  seen  as  the  amount 
of  fadins  increases. 
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Figure  19.  Performance  of  rate  1/2  convolutional  code  with  no  fading:  y2  =  oo  and 
M  =  A  under  optimum  partial-band  jamming  conditions  with  number 
of  hops-pcr-  symbol  (L)  as  a  parameter  when  £4/A"0  =  15  dB. 
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Figure  20.       Performance  of  rate  1/2  convolutional  code  in  Rician  fading:     y2  =  15 

and    M=A   under    optimum   partial-band   jamming    conditions   with 
number  of  hops-per-  symbol  (L)  as  a  parameter  when  EJi\  =  15  dB. 
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Figure  21.      Performance  of  rate  1/2  convolutional  code  in  Rician  fading:     y2  =  5  and 

M  =  4  under  optimum  partial-band  jamming  conditions  with  number 
of  hops-pcr-  symbol  (L)  as  a  parameter  when  EJNn  =  15  dR. 
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Figure  22.      Performance  of  rate  1/2  convolutional  code  in  Rayleigh  fading:     y2  =  0 

and    M  =  A    under    optimum   partial-band   jamming   conditions   with 
number  of  hops-per-symbol  (L)  as  a  parameter  when  EJN0=  15  dB. 
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Figure  23.      Performance  of  rate  1/2  convolutional  code  in  Rician  fading:     y2  =  15 

and    M  =  2   under   optimum   partial-band  jamming  conditions  with 
number  of  hops-pcr-  symbol  (L)  as  a  parameter  when  EJi\0  —  15  dB. 
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Figure  24.      Performance  of  rate  1/2  convolutional  code  in  Rician  fading:     y2  =  5  and 

M  =2  under  optimum  partial-band  jamming  conditions  with  number 
of  hops-per-symbol  (L)  as  a  parameter  when  EJi\0  =  15  dB. 
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Figure  25.       Performance  of  rate  1/3  convolutional  code  in  Rician  fading:     y2  ■■--  15 

and   M  =  4   under   optimum   partial-band   jamming    conditions   with 
number  of  hops-pcr-  symbol  (L)  as  a  parameter  when  EJN0  =  15  dI3. 
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Figure  26.       Performance  of  rate  1/3  convolutional  code  in  Rayleigh  fading:     y2  =  0 

and   M  =  4   under    optimum   partial-band   jamming   conditions   with 
number  ofhops-per-  symbol  (L)  as  a  parameter  when  EJN0  =  15  dB. 
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Figure  27.      Performance  of  Reed-Solomon  (15,8)  code  in  Rician  fading:     y2  =  15  and 

M  =  A  under  optimum  partial-band  jamming  conditions  with  number 
of  hops-pcr-  symbol  (L)  as  a  parameter  when  EbjN0  =  15  dB. 
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Figure  28.       Performance  of  Reed-Solomon  (15,8)  code  in  Rayleigh  fading:     y7  =  0 

and    M  =  4   under    optimum   partial-band   jamming   conditions   with 
number  of  hops-pcr-  symbol  (L)  as  a  parameter  when  EtIN0  =  15  dB. 
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Figure  29.      Performance  comparison  of  codes  with  no  fading:     y7  =  oo  and  .\/  =  4 

under  optimum  partial-band  jamming  conditions  with  L  —  2  (solid  line) 
and  L  =  3  (dashed  line)  when  EJi\  =  15  dB. 
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Figure  30.      Performance  comparison  of  codes  in  Rician  fading:     y2  =  15  and  M  =  4 

under  optimum  partial-band  jamming  conditions  with  L  =  2  (solid  line) 
and  L  =  3  (dashed  line)  when  EJN0  =  15  dB. 
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Figure  31.      Performance  comparision  of  codes  in  Rician  fading:     y2  =  5  and  M  —  A 

under  optimum  partial-band  jamming  conditions  with  L  =  2  (solid  line) 
and  L  =  3  (dashed  line)  when  EJi\0  =  15  dB. 
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Figure  32.       Performance  comparision  of  codes   in   Ravleigh  fading:     y2  =  0  and 

M  =  A  under  optimum  partial-band  jamming  conditions  with   L  =  2 
(solid  line)  and  L  =  3  (dashed  line)  when  £4/A"0  =  15  dB. 
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V.     CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  theoretical  study  and  numerical  simulation  of  randomly  spaced  satellites 
indicate  that  there  will  usually  be  many  satellites  within  communication  range  of 
a  ground  terminal  with  relatively  few  periods  where  coverage  is  not  available. 
The  large  numbers  of  satellites  could  lead  to  high  levels  of  interference  and  many 
multipath  channels  if  no  processing  of  signals  or  directional  transmission  of 
signals  is  performed.  Furthermore,  as  the  range  between  ground  terminals  de- 
siring to  communicate  increases,  the  number  of  satellites  available  to  communi- 
cate through  decreases  exponentially.  This  could  lead  to  longer  periods  of  no 
communication  ability  when  the  satellites  are  not  equipped  with  a  cross-link  ca- 
pability. For  a  250  satellite  array,  from  five  to  20  satellites  can  be  expected 
within  range  about  95%  of  the  time. 

The  theory  needed  to  analyze  the  performance  of  a  fast-frequency-hopped 
MFSK  system  with  side  information  in  a  fading  channel  and  optimum  partial 
band  jamming  has  been  developed.  In  the  configurations  considered,  diversity 
and  coding  both  provide  improvement  in  system  performance.  The  system  is  ex- 
pected to  operate  in  a  fading  channel  and  will  thus  require  diversity  and  coding 
for  improved  performance.  The  optimum  level  of  diversity  for  most  combinations 
of  noise  level  and  fading  Level  is  when  L  =  2  or  L  =  3  hops  per  symbol  are  used. 
In  addition,  the  best  coding  performance  was  obtained  with  the  best  rate  1/2 
convolutional  code.  The  results  for  the  Reed-Solomon  codes  may  not  be  an  ac- 
curate representation  of  their  benefit  in  a  bursty  channel  since  burst  analysis  was 
not  simulated. 

B.  RECOMMENDATIONS 

Further  study  into  other  waveform  configurations  suitable  for  the  LASAT 
system  is  needed.  This  study  could  include  determining  the  benifits  of  direct  se- 
quence spread  spectrum  and  phase  shift  keyed  (PSK)  modulation.   Link  analysis 
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must  also  be  performed  on  the  proposed  LASAT  system.  Multiuser  interference, 
and  a  study  of  antenna  design  should  be  included  in  the  analysis.  The  fast- 
frequency-hopping  MFSK  system  theoretically  developed  in  this  thesis  should 
perform  well  in  a  fading  channel  with  jamming.  To  confirm  the  performance  re- 
sults, a  hardware  simulator  implementing  this  design  should  be  constructed. 

Network  analysis  for  the  LASAT  system  should  also  be  studied  and  should 
consist  of  the  following  tasks: 

•  Allocating  satellite  capacity  to  voice  channels  in  the  case  where  voice  is  not 
packetized. 

•  Design  packet  protocols  to  handle  a  mixture  of   packetized  voice  and  data 

traffic. 

•  Design  highly  reliable  data  protocols  for  priority  messages  to  be  used  when 
the  link  is  forced  to  operate  at  its  threshold. 

•  Design  connectivity  algorithm  for  dvnamic  crosslinks  between  LASAT's. 
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APPENDIX 

The  unconditional  probability  density  function  of  zk  can  be  derived  by  inte- 
grating the  conditional  pdf  over  the  distribution  for  ak  such  that 


fzj-lk)  =         fZlk\Ak(Z^  '  "kVAftkWk- 


(AA) 


where  fz  \A{zik\ak)  and  fAk{ak)  are  given  by  (4.17)  and  (4.4),  respectively.   There- 
fore, the  integral  to  be  solved  is 
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The  solution  of  the  integral  may  be  found  by  using  the  following  indentity  from 
Ref.  17:  p.  314 
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where   Jv(v)   is   the    Bessel    function    of   the    first    kind.      If  it    is   noted    that 
I0(x)  =  J0(/A").  and  following  substitutions  are  made 
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the  result  is 
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This  equation  can  be  put  into  a  more  convienent  format  if  the  substitutions  of 
ph  =  ol\R  =  aljlal  and  /?  =  2cr2^  =  2g2\2u\  are  made 
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where  /),.  is  the  output  sienal-to-noise  ratio  of  the  specular  (direct)  component, 
and  ft  is  the  output  signal-to-  noise  ratio  of  the  diffuse  component. 
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